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Abstract 


The unique ultra-relativistic, massless, nature of electron states in two-dimensional ex¬ 
tended graphene sheets, brought about by the honeycomb lattice arrangement of carbon atoms 
in two-dimensions, provides ingress to explorations of fundamental physical phenomena in 
graphene nanostructures. Here we explore the emergence of new behavior of electrons in 
atomically precise segmented graphene nanoribbons (GNRs) and graphene rings with the use 
of tight-binding calculations, non-equilibrium Green’s function transport theory, and a newly 
developed Dirac continuum model that absorbs the valence-to-conductance energy gaps as 
position-dependent masses, including topological-in-origin mass-barriers at the contacts be¬ 
tween segments. Through transport investigations in variable-width segmented GNRs with 
armchair, zigzag, and mixed edge terminations we uncover development of new Fabry-Perot- 
like interference patterns in segmented GNRs, a crossover from the ultra-relativistic mass¬ 
less regime, characteristic of extended graphene systems, to a massive relativistic behavior 
in narrow armchair GNRs, and the emergence of nonrelativistic behavior in zigzag-terminated 
GNRs. Evaluation of the electronic states in a polygonal graphene nanoring under the influence 
of an applied magnetic field in the Aharonov-Bohm regime, and their analysis with the use of 
a relativistic quantum-field theoretical model, unveils development of a topological-in-origin 
zero-energy soliton state and charge fractionization. These results provide a unifying frame¬ 
work for analysis of electronic states, coherent transport phenomena, and the interpretation of 
forthcoming experiments in segmented graphene nanoribbons and polygonal rings. 


Introduction 

In the last three deeades transport through moleeular junetions^^^has attraeted mueh attention be- 
eause of fundamental aspeets of the proeesses involved, as well as of potential praetieal prospeets. 
In partieular, studies in this direetion have intensified sinee the diseovery of new forms of earbon 
allotropes, starting with the fullerenes^^ and earbon nanotubes^ (CNTs) in the 1980s and 1990s 
respeetively, and the isolation of graphene^^ in 2004. The above earbon allotropes differ in shape 
(eurvature), topology, and dimensionality, with the fullerenes being zero-dimensional (OD) with 
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spherical or prolate shape, the CNTs being one-dimensional (ID) cylinders, and graphene be¬ 
ing a two-dimensional (2D) plane (or ID planar ribbons^^^®^. In the fullerenes the carbon-atom 
network is made of (non-adjacent) hexagons and pentagons, whereas the CNTs and graphene are 
entirely hexagonal lattices (curved in CNTs) described in terms of a unit cell with a two-atom basis 
with the two carbon atoms occupying two sublattices (A and B, also mapped into the up and down 
pseudospin states^. In the absence of defects, in-plane (c) bonding occurs through sp^ hybrid 
orbitals and out-of-plane bonding {n) involves the orbital; in the following only the physics of 
;r-states is considered. 

The hexagonal network topology of graphene gives rise to relativistic behavior of the low- 
energy excitations which is captured by the ultra-relativistic massless Dirac-Weyl (DW) equation, 
with the Fermi velocity of graphene {vp = c/300) replacing the velocity of light. ^ Among the 
many manifestations of the relativistic behavior in graphene is the Klein paradox, that is "... unim¬ 
peded penetration of relativistic particles through high and wide potential barriers - is one of the 
most exotic and counterintuitive consequences of quantum electrodynamics.”^ The surprising rel¬ 
ativistic behavior in graphene has indeed been recognized in the 2010 Nobel award in physics to 
A. Geim and K. Novosolev. 

Another carbon-based system that was the subject of an earlier (2000) chemistry Nobel award 
(to A.J. Heeger, A.G. MacDiarmid, and H. Shirakawa) is conductive polymers, with polyacetelyne 
(PAG) being a representative example.^ PAG is a ID chain of carbon atoms forming a conjugated 
polymer with sp^p^ hybridization that leads to one unpaired electron per carbon atom (half-filled n- 
band) as in graphene. Linearizing the spectrum of this ID equally-spaced carbon chain at the Fermi 
level (that is at the Dirac points, i.e., the zero-energy points of the energy vs momentum dispersion 
relation) results in a ID massless Dirac-equation description of the low energy excitations of the 
system.fi^E^ 

However, the equally-spaced ID system is unstable and consequently it distorts (structurally) 
spontaneously (Peierls distortion!^ yielding a modulation (alternation) of the spacing between 
successive sites that results in dimerization of successive atoms along the chain and the opening of 


3 




a gap in the electronic spectrum. This dimerization can occur in two energetically degenerate, but 
spatially distinct, patterns, termed as equivalent domain structures. Either of these domains is a 
realization of the ID Dirac equation with a constant mass term, ^. Disruption of the dimerization 
pattern (e.g., by transforming from one domain to another along the chain) creates a domain wall 
which can be described in the ID Dirac formulation through the use of a position-dependent mass 
term of alternating sign (with the spatial mass-profile connecting with The solution of 

this generalized ID Dirac equation is a soliton characterized by having zero energy and by being 
localized at the domain wall. In this paper, the Dirac equation with position-dependent mass will be 
used in investigations of the electronic structure, transport, and magnetic-field-induced phenomena 
(Aharonov-Bohm) in graphene nanostructures such as nanoribbon junctions and rings.^^^^ 

From the above we conclude that in the two extreme size-domains, that is a 2D infinite graphene 
sheet and a ID carbon chain (PAC), the systems exhibit behavior that is relativistic in nature. In an 
attempt to bridge between these two size-domains, we briefly discuss in the following systems of 
successively larger width, starting from the polyacene (a quasi-ID chain of fused benzene rings) 
which may be regarded as the narrowest graphene nanoribbon (GNR) with zigzag edge termina¬ 
tions. 

Polyacene was investigated^ first in 1983. It was found that, as in graphene and PAC (without 
distortion), the valence and conduction bands of undistorted polyacene touch at the edge of the 
Brillouin zone. However, unlike 2D graphene and PAC, the dispersion relation about the touching 
point is quadratic, conferring a non-relativistic (Schrodinger equation) character. We show in this 
paper that this surprising finding persists for sufficiently narrow GNRs with zigzag termination 
(zGNRs). However, narrow armchair-terminated GNRs (aGNRs) are found here to maintain rela¬ 
tivistic behavior, with metallic ones being massless and semiconducting ones being massive (both 
classes obeying the Einstein energy relation). 

It this paper, we discuss mainly manifestations of relativistic and/or nonrelativistic quantum 
behavior explored through theoretical considerations of transport measurements in segmented 
graphene nanoribbons of variable width, and spectral and topological effects in graphene nanorings 
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in the presence of magnetic fields. 

Transport through narrow graphene channels — particularly bottom-up fabricated and atomically- 
precise graphene nanoribbons^^^^ — is expected to offer ingress to unique behavior of Dirac 
electrons in graphene nanostructures. In particular, the wave nature of elementary particles (e.g., 
electrons and photons) is commonly manifested and demonstrated in transport processes. Because 
of an exceptionally high electron mobility and a long mean-free path, it has been anticipated that 
graphene)^ devices hold the promise for the realization, measurement, and possible utilization 
of fundamental aspects of coherent and ballistic transport behavior, which to date have been ob¬ 
served, with varying degrees of success, mainly at semiconductor interfacesquantum point 
contacts,® metallic wires,® and carbon nanotubes.® 

Another manifestation of coherent ballistic transport are interference phenomena, reflecting 
the wave nature of the transporting physical object, and associated most often with optical (elec¬ 
tromagnetic waves, photons) systems or with analogies to such systems (that is, the behavior of 
massless particles, as in graphene sheets). Measurements of interference patterns are commonly 
made with the use of interferometers, most familiar among them the multi-pass optical Fabry-Perot 
(OFP) interferometer.®® The advent of 2D forms of carbon allotropes has motivated the study of 
optical-like interference phenomena associated with relativistic massless electrons, as in the case 
of metallic carbon nanotubes® and graphene 2D p-n junctions.® (We note that the hallmark of 
the OFP is that the energy separation between successive maxima of the interference pattern varies 
as the inverse of the cavity length L.) 

For GNRs with segments of different widths, our investigations reveal diverse Fabry-Perot 
transport modes beyond the OFP case, with conductance quantization steps {nGo,n = 1,2,3,..., 
with Go = jh) found only for uniform GNRs. In particular, three distinct categories of Fabry- 
Perot interference patterns are identified: 

1. FP-A: An optical FP pattern corresponding to massless graphene electrons exhibiting equal 
spacing between neighboring peaks. This pattern is associated with metallic armchair nanorib¬ 
bon central segments. This category is subdivided further to FP-Al and FP-A2 depending 
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on whether a valenee-to-eonduetion gap is absent (FP-Al, assoeiated with metallie armehair 
leads), or present {FP-A2, eorresponding to semieondueting armehair leads). 

2. FP-B: A massive relativistic FP pattern exhibiting a shift in the eonduetion onset due to 
the valenee-to-eonduetion gap and unequal peak spacings. This pattern is assoeiated with 
semieondueting armehair nanoribbon eentral segments, irrespeetive of whether the leads are 
metallie armehair, semieondueting armehair, or zigzag. 

3. FP-C: A massive non-relativistic FP pattern with 1/L^ peak spacings, but with a vanishing 
valenee-to-eonduetion gap, L being the length of the eentral segment. This pattern is the one 
expeeted for usual semieonduetors deseribed by the (nonrelativistie) Sehrodinger equation, 
and it is assoeiated with zigzag nanoribbon eentral segments, irrespeetive of whether zigzag 
or metallie armehair leads are used. 

We report in this paper on the unique apeets of transport through segmented GNRs obtained 
from tight-binding non-equilibrium Green’s funetionl^^lUJ (TB-NEGF) ealeulations in eonjunetion 
with an analysis based on a one-dimensional (ID) relativistie Dirae eontinuum model. This eontin- 
uum model goes beyond the physies of the massless Dirae-Weyl (DW) eleetron, familiar from two- 
dimensional (2D) honeyeomb earbon sheets,^ and it is referred to by us as the Dirae-Fabry-Perot 
(DFP) theory (see below for the ehoiee of name). In partieular, it is shown that the valenee-to- 
eonduetion energy gap in armehair GNR (aGNR) segments, as well as the barriers at the interfaees 
between nanoribbon segments, ean be ineorporated in an effeetive position-dependent mass term 
in the Dirae hamiltonian; the transport solutions assoeiated with this hamiltonian exhibit eondue- 
tanee patterns eomparable to those obtained from the mieroseopie NEGF ealeulations. For zigzag 
graphene nanoribbon (zGNR) segments, the valenee-to-eonduetion energy gap vanishes, and the 
mass term is eonsonant with exeitations eorresponding to massive nonrelativistie Sehrodinger-type 
earriers. The faithful reproduetion of these unique TB-NEGF eonduetanee patterns by the DFP the¬ 
ory, ineluding mixed armehair-zigzag eonfigurations (where the earriers transit from a relativistie 
to a nonrelativistie regime), provides a unifying framework for analysis of eoherent transport phe- 
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nomena and for interpretation of experiments targeting fundamental understanding of transport in 
GNRs and the future development of graphene nanoelectronics. 

To demonstrate the aforementioned soliton formation due to structural topological effects (dis¬ 
cussed by us above in the context of poly acetylene), we explore with numerical tight-binding 
calculations and a Dirac-Kronig-Penney (DKP) approach, soliton formation and charge fraction- 
ization in graphene rhombic rings; this DKP approach is based on a generalized Dirac equation 
with alternating-sign position-dependent masses. 

Before leaving the Introduction, we mention that, due to their importance as fundamental phe¬ 
nomena, Aharonov-Bohm-type effects in graphene-nanoribbons systems have attracted (in ad¬ 
dition to Refs. I20l - l2^ considerable theoretical attention. These latter theoretical papers, 
however, based their analysis exclusively on tight-binding and/or DFT calculations,l®^l^®^I^ 
or they used in addition a two-dimensional Dirac equation with infinite-mass boundary condi- 
tions .03^ Transcending the level of current understanding which explores direct similaritiess with 
the Aharonov-Bohm physics in semiconductor and metallic mesoscopic rings, our work here ana¬ 
lyzes the TB results in conjunction with a continuum ID generalized Dirac equation (that incorpo¬ 
rates a position-dependent mass term), and thus it enables investigations of until-now unexplored 
topological aspects and relativistic quantum-field analogies of the AB effect in graphene nanosys¬ 
tems. 

Furthermore we note that oscillations in the conductance of graphene nanoribbons in the pres¬ 
ence of magnetic barriers were found in a theoretical studyl^ (using exclusively a TB-NEGF ap¬ 
proach), as well as in an experimental investigation of high quality bilayer nanoribbons,® and they 
were attributed to Fabry-Perot-type interference. In the absence of a continuum Dirac analysis, the 
precise relation of such oscillatory patterns to our Fabry-Perot theory (based on the incorporation 
of a position-dependent mass term in the Dirac equation at zero magnetic field) warrants further 
investigation. 

Finally, gap engineering in graphene ribbons under strong external fields was studied in Ref. 
[501 We stress again that one of the main results in this paper is the appearance of “forbidden” 
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solitonic states inside the energy gap in the context of the low-magnetic-field Aharonov-Bohm 
spectra of graphene nanorings; see the part titled “Aharonov-Bohm spectra of rhombic graphene 
rings” in the Results section. 

Methods 

Dirac-Fabry-Perot model 

The energy of a fermion (with one-dimensional momentum px) is given by the Einstein relativis¬ 
tic relation E = ^ {PxVfY + {-y^VpY, where ^ is the rest mass and is the Fermi velocity of 
graphene. In a gapped uniform graphene nanoribbon, the mass parameter is related to the particle- 
hole energy gap, A, as = ^/(2v|.). Following the relativisitic quantum-field Fagrangian formal¬ 
ism, the mass ^ is replaced by a position-dependent scalar Higgs field ^{x) = m{x)vp, to which 
the relativistic fermionic field 'F(x) couples through the Yukawa Fagrangian^ .ify = — 0'F^/3'F 
(j8 being a Pauli matrix). For ^{x) = (constant) ^Vp = ^o, and the massive fermion Dirac 
theory is recovered. The Dirac equation is generalized as (here we do not consider applied electric 
or magnetic fields) 

[E-V{x)]^' + ihvfa-^-l5^{x)^^ = 0. (1) 

In one dimension, the fermion field is a two-component spinor 'P = (t/A„, t//;)^; u and I stand, 
respectively, for the upper and lower component and a and /3 can be any two of the three Pauli 
matrices. Note that the Higgs field enters in the last term of Eq. (??). V{x) in the first term is 
the usual electrostatic potential, which is inoperative due to the Klein phenomenon and thus is 
set to zero for the case of the armchair nanoribbons (where the excitations are relativistic). The 
generalized Dirac Eq. (??) is used in the construction of the transfer matrices of the Dirac-Fabry- 
Perot model described below. 

The building block of the DFP model is a 2 x 2 wave-function matrix ■ formed by the compo¬ 
nents of two independent spinor solutions (at a point x) of the onedimensional, first-order gener- 
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alized Dirac equation [see Eq. (3) in the main paper]. ■ plays^the role of the Wronskian matrix 
W used in the seeond-order nonrelativistie Kronig-Penney model. Following Ref. [5T] and gener¬ 
alizing to N regions, we use the simple form of ■ in the Dirae representation (a = Ci, /3 = ( 73 ), 
namely 

^ ^iKx ^—iKx 


■x(^) = 


( 2 ) 


V 


j^eiKx 




/ 


where 




2,,4 


m V 


F 




A = 


fivpK 


E — V + mvp 


(3) 


The transfer matrix for a given region (extending between two matehing points xi and X 2 speeifying 
the potential steps m-”^) is the produet M/i:(ji;i,X 2 ) = ■x( 2 C 2 )*^^ (-^ 1 ); this latter matrix depends only 
on the width X 2 —xi of the region, and not separately on xi or X 2 . 

The transfer matrix eorresponding to a series of N regions ean be formed^ as the produet 


(4) 

i=l,N 

where — x,| = L,- is the width of the ith region [with (m,y,^,A) —)■ (m,, The 

transfer matrix assoeiated with the transmission of a free fermion of mass ^ (ineoming from the 
right) through the multiple mass barriers is the produet 


^{E) =m,^\xN+i)ti,N+imk{xi), (5) 

W\i\ik= ^{E-VY-J^^v‘],/{hvF), \E-V\ > \ for armehair leads V = 0, while for zigzag 

leads V = Naturally, in the ease of metallie armehair leads, k = E/(hvp), sinee ^ = 0. 

Then the transmission eoeffieient T is 


^ l^22P’ 


( 6 ) 
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while the refleetion eoeffieient is given by 


R = 



(7) 


At zero temperature, the eonduetanee is given by G = {2e^/h)T ; T is the transmission eoeffi- 
eient in Eq. (??). 


Dirac-Kronig-Penney superlatice model 

The transfer matrix eorresponding to either half of the rhombus ean be formed^^ as the produet 

t„ = ]^ xi = 0, JC4 = L, (8) 

(= 1,3 

with L being the length of half of the perimeter of the rhombus; L = Li +L 2 + L 3 , with Li = L 3 = a 
and L 2 = b. The transfer matrix assoeiated with the eomplete unit eell (eneireling the rhombie 
ring) is the produet 

2 

T=nt- (9) 

n=l 

Following Refs. 1211521 we eonsider the superlattiee generated from the virtual periodie trans¬ 
lation of the unit eell as a result of the applieation of a magnetie field B perpendieular to the ring. 
Then the Aharonov-Bohm energy speetra are given as solutions of the dispersion relation 

eos [2;r(4>/d>o + 77 )]= Tr[T(E)]/2, (10) 

where we have explieitly denoted the dependenee of the r.h.s. on the energy E; rj = 0 for a rhombus 
with type-1 eorners and 7] = 1/2 for a rhombus with type-11 eomers. 

The energy speetra and single-partiele densities do not depend on a speeifie representation. 
However, the wave funetions (upper and lower spinor eomponents of the fermionie field 'E) do 
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depend on the representation used. To transform the initial DKP wave functions to the (a = 02 , 
P = Oi) representation, which corresponds to the natural separation of the tight-binding amplitudes 
into the A and B sublattices, we apply successively the unitary transformations D23 = {(72 + Cs) / a /2 
and D3 — exp{ina3 / 4). 

TB-NEGF formalism 

To describe the properties of graphene nanostructures in the tight-binding approximation, we use 
the hamiltonian 

HTB = -t cjcj + h.c., (11) 

<i,j> 

with <> indicating summation over the nearest-neighbor sites ij. t = 2.7 eV is the hopping 
parameter of two-dimensional graphene. 

To calculate the TB-NEGF transmission coefficients, the Hamiltonian (??) is employed in con¬ 
junction with the well known transport formalism which is based on the nonequilibrium Green’s 
functions 

According to the Landauer theory, the linear conductance is G{E) = {2e^/h)T{E), where the 
transmission coefficient is calculated as T{E) = Tx\r The Green’s function ‘^{E) is 
given by 

<^{E) = (E + ir] - //d- _^ (12) 

with //'pb' being the Hamiltonian of the isolated device (junction without the leads). The self¬ 
energies are given by where the hopping matrices Tp(^) 

describe the left (right) device-to-lead coupling, and ^ is the Hamiltonian of the semi-infinite 
left (right) lead. The broadening matrices are given by 
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Figure 1: Conductance quantization steps (a) for a uniform semiconducting armchair 
nanoribbon (I) contrasted to Fabry-Perot oscillations (b-f) of two 3-segment armchair GNRs 
[(II) and (IV)] with both a semiconducting central constriction and semiconducting leads 
(I3-7-I3 and 7-13-7). For the systems shown here all the segments have armchair edge termina¬ 
tion (hence, AAA), and all have width corresponding to semiconducting GNRs [hence (sss)]. (Ill, 
V) Schematics of the mass barriers used in the DFP modeling, with the dashed line denoting the 
zero mass. The physics underlying such a junction is that of a massive relativistic Dirac fermion 
impinging upon the junction and performing multiple reflections (above m\Vp) within a particle 
box defined by the double-mass barrier. (c,e) TB-NEGF conductance as a function of the Fermi 
energy of the massive Dirac electrons in the leads, (d) DFP conductance reproducing [in the energy 
range of the IGq step, see (b)] the TB-NEGF result in (c). The mass-barrier parameters used in 
the DFP reproduction were Li = 55ao, miv^ = Q.22t, L 2 = lag, m 2 V^ = 0.5t. The mass of the 
electrons in the leads was ^iVp = 0.166t. (f) DFP conductance reproducing the TB-NEGF result 
in (e). The parameters used in the DFP reproduction were Li = 53.bag, miVp = 0.166t, L 2 = lao, 
m 2 Vp = 0.5It. The mass of the electrons in the leads was ^iVp = 0.22t. gq = 0.246 nm is the 
graphene lattice constant; t = 2.7 eV is the hopping parameter. 
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Results 


Segmented Armchair GNRs: All-semiconducting 

Our results for a 3-segment , where is the lead width and is the 

width of the eentral segment) all-semiconducting aGNR, AAA (sss), are portrayed in Fig. [see 
sehematie lattiee diagrams in Figs.[2l) and[2ll)]. A uniform semieondueting armehair GNR [see 
Fig.[2l)] exhibits ballistie quantized-eonduetanee steps [see Fig. [TJa)]. In eontrast, eonduetanee 
quantization is absent for a nonuniform 3-segment (13 — 7 — 13) aGNR; see Figs, [^b) — [^d)- 
Here oseillations appear instead of quantized steps. The first oseillation appears at an energy 
~ O.llt [Fig. [^b)], whieh refleets the intrinsie gap A/2 of the semieondueting eentral segment 

belonging to the elass II of aGNRs, speeified^i^^i^^^ by a width = 3/ -1- I, / = 1,2,3,_We 

reeall that as a funetion of their width, , the armehair graphene nanoribbons fall into three 
elasses: (I) = 3/ (semieondueting, A > 0), (II) = 3/3-1 (semieondueting, A > 0), and 

(III) = 3/ + 2 (metallie, A = 0), / = 1,2,3,.... 

That the leads are semieondueting does not have any major effeet. This is due to the faet that 
, and as a result the energy gap m\v\ of the eentral segment is larger than the energy 
gap of the semieondueting leads [see sehematie in Fig.[^III)]. In the opposite ease (eentral 
segment wider than the leads), the energy gap of the semieondueting leads would have determined 
the onset of the eonduetanee oseillations. 

The armehair GNR ease with interehanged widths (i.e., 7 — 13 — 7 instead of 13 — 7 — 13) is 
portrayed in Figs.[2e) — [^f)- this ease the energy gap of the semieondueting leads (being the 
largest) determines the onset of the eonduetanee oseillations. It is a testimonial of the eonsisteney 
of our DFP method that it ean reproduee [see Fig. [^d) and Fig. [^f)] both the 13 — 7—13 and 
7 — 13 — 7 TB-NEGF eonduetanees; this is aehieved with very similar sets of parameters taking 
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into consideration the eentral-segment-leads interehange. We note that the larger spaeing between 
peaks (and also the smaller number of peaks) in the 7 — 13 — 7 ease is due to the smaller mass of 
the eentral segment (0.166t instead of 0.22t). 

Further insight ean be gained by an analysis of the diserete energies assoeiated with the humps 
of the eonduetanee oseillations in Fig.[2e) and the resonant spikes in Fig.[2e). Indeed a simplified 
approximation for the eleetron eonfinement in the eontinuum model eonsists in eonsidering the 
graphene eleetrons as being trapped within a ID infinite-mass square well (IMSW) of length Li 
(the mass terms are infinite outside the interval Li and the eoupling to the leads vanishes). The 
diserete speetrum of the eleetrons in this ease is given^^by 

En = + (13) 

where the wave numbers are solutions of the transeendental equation 

tan(k„Li) = —). (14) 


When ^ = 0 [massless Dirae-Weyl (DW) eleetrons], one finds for the speetrum of the IMSW 
model: 

En = {n+l/2)TZflVF/Lu ( 15 ) 

with n = 0,1,2,.... 

E„-En^i=2Eo, n= 1,2,3,..., (16) 

whieh is twiee the energy 

Eo = nhvF/{2Li) (17) 

of the lowest state. 

As is well known, a eonstant energy separation of the intensity peaks, inversely proportional to 
the length of the resonating eavity [here Li, see Eqs. (??) and (??) above] is the hallmark of the 
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optical Fabry-Perot, reflecting the linear energy dispersion of the photon in optiesl^ or a massless 
DW eleetron in graphene struetures.l^Most revealing is the energy offset away form zero of the 
first eonduetanee peak, whieh equals exaetly one-half of the eonstant energy separation between 
the peaks. In one-dimension, this is the hallmark of a massless fermion subjeet to an infinite- 
mass-barrier eonfinement.l^ Naturally, in the ease of a semieondueting segment (see below), this 
equidistant behavior and 1/2-offset of the eonduetanee peaks do not apply; this ease is aeeounted 
for by the Dirae-Fabry-Perot model presented in Methods, and it is more general than the optieal 
Fabry-Perot theory assoeiated with a photonie eavity.^^ 

In the nonrelativistie limit, i.e., when hkn <C ^Vf, one gets 


tan(k„Li) ~ 0, 


(18) 


whieh yields the well known relations mt and 

En ~ ^vl + T?h^TZ^/{2Jii:L\). (19) 

For a massive relativistie eleetron, as is the ease with the semieondueting aGNRs in this paper, one 
has to numerieally solve Eq. (??) and then substitute the eorresponding value of in Einstein’s 
energy relation given by Eq. (??). 

Erom an inspeetion of Eig.[^ one ean eonelude that the physies assoeiated with the all-semieondueting 
AAA junetion is that of multiple refleetions of a massive relativistie Dirae fermion bouneing baek 
and forth from the edges of a “quantum box” ereated by a double-mass barrier [see the sehematie of 
the double-mass barrier in Eigs.[2lII) and[2V)] . In partieular, to a good approximation the energies 
of the eonduetanee oseillation peaks are given by the IMSW Eq. (??) with = miv^ = 0.22t 
(13 — 7 — 13) or = mivjr = 0.166t (7 — 13 — 7). In this respeet, the separation energy be¬ 
tween sueeessive peaks in Eigs. [^b), [2^)’ [3^)’ ^iid[2f) is not a eonstant, unlike the ease of an 
all-metallie junetionl^ (or a photonie eavity^. 

The eonduetanee patterns in Eigs. [^d) and[2f) eorrespond to the eategory FP-B (see the in- 
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troductory section). These patterns cannot be accounted for by the optical Fabry-Perot theory, 
but they are well reproduced by the generalized Dirac-Fabry-Perot model introduced by us in the 
Methods section. 
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Figure 2: (Left) Conductance for a 3-segment nanoribbon, all segments having armchair edge 
termination (i.e. AAA), with a metallic = 5) central constriction and semiconducting 
leads = 7), hence the designation (sss); see schematic lattice diagram in (I). (II) Schematics 
of the mass barriers used in the DFP modeling, (a) TB-NEGF conductance as a function of the 
Fermi energy of the massive Dirac electrons in the leads, (b) DFP conductance reproducing the TB- 
NEGF result in (a). The mass-barrier parameters used in the DFP reproduction were Li = 60.4ao, 
mi = 0, L 2 = lao, m 2 Vp = 031t. The mass of the electrons in the leads was = 0.23t. (c)-(d) 
The total DOS of the junction and the density of states in the isolated leads, respectively, according 
to the TB-NEGF calculations. The arrows indicate the onset of the electronic bands in the leads. 
Note that the DOS in (c) reveal the existence of five sharp electronic states below the onset (at 
0.23t = ^iv^) of the first band in the leads [see (d)], which consequently do not generate any 
conductance resonances [see (a) and (b)]. Note further in (c) the equal energy spacing between the 
vertical lines [the five solid (red) and four dashed (black) ones] associated with the resonances of 
a massless electron confined within the central metallic aGNR segment. 

(Right) H-passivation effects in the conductance of a 3-segment armchair nanoribbon with 
a metallic (^'^ = 5) central constriction and semiconducting leads {3^^ = 7); see schematic 
lattice diagram in (III). Note that the nearest-neighbor C-C bonds at the armchair edges (thick red 
and blue lines) have hopping parameters t' = 1.12t. (IV) Schematics of the position-dependent 
mass field used in the DFP modeling, (e) TB-NEGF conductance as a function of the Fermi 
energy, (f) DFP conductance reproducing the TB-NEGF result in (e). The mass parameters used 
in the DFP reproduction were Li = 59.5ao, m^v^ = O.OSt, L 2 = l.Sao, m 2 v\ = 0.30t. The mass of 
the electrons in the leads was = 0.28t. (g)-(h) The total DOS of the junction and the density 
of states in the isolated leads, respectively, according to the TB-NEGF calculations. The arrows 
indicate the onset of the electronic bands in the leads; note the shifts from 0.23t to 0.28t and from 
0.42t to 0.38t for the onsets of the first and second bands, respectively, compared to the case with 
t' = t m (d). Compared to left part of Fig.[^ the subtle modifications of mass parameters brought 
about by having t' = 1.12t result in having six sharp electronic states [see (g)] below the onset 
(at 0.28t = ^iv^) of the first band in the leads [see (h)], which consequently do not generate any 
conductance resonances [see (e) and (f)]. In addition, within the energy range (0 to 0.40 shown 
in (e) and (f), there are now only two conducting resonances, instead of three compared to (a) and 
(b). ao = 0.246 nm is the graphene lattice constant; t = 2.7 eV is the graphene hopping parameter. 
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Segmented Armchair GNRs: semiconducting-metallic-semiconducting. 

Our results for a 3-segment (7 — 5 — 7) semiconducting-metallic-semiconducting aGNR, AAA 
(sms), are portrayed in Fig. [left, see sehematie lattiee diagram in Fig. |^I)]. The first FP os- 
eillation in the TB-NEGF eonduetanee displayed in Fig.|^a) appears at an energy ~ 0.22t, whieh 
refleets the intrinsie gap A/2 of the semieondueting leads (with = 7). The energy spaeing 
between the peaks in Fig. [^a) is eonstant in agreement with the metallie (massless DW eleetrons) 
eharaeter of the eentral segment with = 5. The TB-NEGF pattern in Fig.|^a) eorresponds to 
the Fabry-Perot eategory FP-A2. As seen from Fig.|^b), our generalized Dirae-Fabry-Perot theory 
is again eapable of faithfully reprodueing this behavior. 

A deeper understanding of the AAA (sms) ease ean be gained via an inspeetion of the density 
of states (DOS) plotted in Fig. |^e) for the total segmented aGNR (eentral segment plus leads) 
and in Fig. [^d) for the the isolated leads. In Fig. [^e), nine equidistant resonanee lines are seen. 
Their energies are elose to those resulting from the IMSW Eq. (??) (with Li = 60.4ao, see the 
eaption of Eig. for a massless DW eleetron. Out of these nine resonanees, the first five do not 
eonduet [eompare Eigs. |^a) andj^e)] beeause their energies are lower than the minimum energy 
(i.e., A/2 = ^iVp ~ 0.230 of the ineoming eleetrons in the leads [see the onset of the first band 
(marked by an arrow) in the DOS eurve displayed in Eig.j^d)]. 

Segmented Armchair GNRs: Effects of hydrogen passivation. 

As shown in Refs. \nm a detailed deseription of hygrogen passivation requires that the hop¬ 
ping parameters t' for the nearest-neighbor C-C bonds at the armehair edges be given by t' = 
1.12t. Taking this modifieation into aeeount, our results for a 3-segment semiconducting-metallic- 
semiconducting aGNR are portrayed in Eig. [right, see sehematie lattiee diagram in Eig.|^III)]; 
this lattiee eonfiguration is denoted as "AAA (sms) H-passivation." The first EP oseillation in the 
TB-NEGE eonduetanee displayed in Eig. [^e) appears at an energy ~ 0.28t, whieh refleets the 
intrinsie gap A/2 of the properly passivated semieondueting leads (with = 7). The energy 
spaeing between the peaks in Eig. |^e) is slightly away from being eonstant in agreement with the 
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small mass miVp = 0.05f acquired by the central segment with = 5, due to taking t' = 1.12t. 
As seen from Fig. |^f), our generalized Dirac-Fabry-Perot theory is again capable of faithfully 
reproducing this behavior. 

A deeper understanding of the AAA (sms)-H-passivation case can be gained via an inspection 
of the DOS plotted in Fig. |^g) for the total segmented aGNR (central segment plus leads) and in 
Fig.j^h) for the isolated leads. In Fig.j^g), eight (almost, but not exactly, equidistant) resonance 
lines are seen. Their energies are close to those resulting from the IMSW Eq. (??) (with L\ = 
59.5ao and miVp = 0.05t; see the caption of Fig. for a Dirac electron with a small mass. Out 
of these eight resonances, the first six do not conduct [compare Figs, [^e) and|^g)] because their 
energies are lower than the minimum energy (i.e., A/2 = ^iVp ~ 0.280 of the incoming electrons 
in the leads [see the onset of the first band (marked by an arrow) in the DOS curve displayed in 
Fig.|^h)]. From the above we conclude that hydrogen passivation of the aGNR resulted in a small 
shift of the location of the states, and opening of a small gap for the central metallic narrower 
(with a width of = 5) segment, but did not modify the conductance record in any qualitative 
way. Moreover, the passivation effect can be faithfully captured by the Dirac FP model by a small 
readjustment of the model parameters. 

All-zigzag segmented GNRs. 

It is interesting to investigate the sensitivity of the interference features on the edge morphology. 
We show in this section that the relativistic transport treatment applied to segmented armchaie 
GNRs does not maintain for the case of a nanoribbon segment with zigzag edge terminations. In 
fact zigzag GNR (zGNR) segments exhibit properties akin to the well-known transport in usual 
semiconductors, i.e., their excitations are governed by the nonrelativistic Schrodinger equation. 

Before discussing segmented GNRs with zigzag edge terminations, we remark that such GNRs 
with uniform width exhibit stepwise quantization of the conductance, similar to the case of a 
uniform armchair-edge-terminated GNR [see Fig.[2a)]. 

In Fig. I^a), we display the conductance in a three-segment junction [see lattice schematic in 
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Figure 3: Conductance for ZZZ (all-zigzag edge termination, left column) and AZA 
(armchair-zigzag-armchair edge termination, right column) segmented nanorihhon junc¬ 
tions. See corresponding lattice diagrams in (I) and (II). The 3 -segment GNRs are denoted as 
, with (i =1,2) being the number of carbon atoms specifying the width 
of the ribbon segments. The armchair leads in the AZA junction are metallic = 23, class 
III aGNR). (a)-(b) TB-NEGF conductance for the ZZZ and AZA junction, respectively, (c) DFP 
conductance reproducing the TB-NEGF result in (a) for the ZZZ junction, (d) DFP conductance 
reproducing the TB-NEGF result in (b) for the AZA junction. In spite of the different edge mor¬ 
phology, the Fabry-Perot patterns in (a) and (b) are very similar. The central zigzag segment con¬ 
trols the Fabry-Perot patterns. According to the continuum DFP analysis, the physics underlying 
such patterns is that of a massive nonrelativistic Schrodinger fermionic carrier performing multiple 
reflections within a cavity defined by a double-mass barrier [see diagram in (III)], but with the ad¬ 
ditional feature that Vi = —m\Vp and V/ = —^iVp are also considered for segments or leads with 
zigzag edge terminations (see text for details). The mass and V, parameters used in the DFP calcula¬ 
tions were Li = 30ao, rniv\ = 2.231 — cEt, with c = 7.3, Vi = —miv^, L 2 = l.lao, m 2 Vp = 0.38t, 
V 2 = —m2Vp/3, .y£iVp = 2.30t, Vj = —.y^iVp in (c) and Li = 29.1ao, = 2.65t — cEt, with 

c = 8.4, Vi = —miVp, L 2 = l.Oao, m 2 Vp = 0.30t, V 2 = —m 2 Vp, .Mi = 0, V/ = 0 in (d). Mi and V/ 
denote parameters of the leads. E is the energy in units oft. ao = 0.246 nm is the graphene lattice 
constant; t = 2.7 eV is the hopping parameter. 
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Fig.|^I)] when all three segments have zigzag edge terminations (denoted as ZZZ), but the eentral 
one is narrower than the lead segments. The main finding is that the eentral segment behaves again 
as a resonant eavity that yields an oseillatory eonduetanee pattern where the peak spaeings are un¬ 
equal [Fig.j^a)]. This feature, whieh deviates from the optieal Fabry-Perot behavior, appeared also 
in the DFP patterns for a three-segment armehair junetion whose eentral segment was semieonduet- 
ing, albeit with a different dependenee on L [see Figs. and Eq. (??)]. Moreover, from a set of 
systematie ealeulations (not shown) employing different lengths and widths, we found that the en¬ 
ergy of the resonant levels in zGNR segments varies on the average as ~ (n/L)^, where the integer 
n eounts the resonanees and L indieates the length of the eentral segment. However, a determining 
differenee with the armehair GNR ease in Fig. is the vanishing of the valenee-to-eonduetanee 
gap in the zigzag ease of Fig. |^a). It is well known that the above features are assoeiated with 
resonant transport of eleetronie exeitations that obey the nonrelativistie seeond-order Sehrodinger 
equation. 

Naturally, one eould formulate a eontinuum transport theory based on transfer matriees (see 
Methods) that use the ID Sehrodinger equation instead of the generalized Dirae Eq. (??). Sueh a 
Sehrodinger-equation eontinuum approaeh, however, is unable to deseribe mixed armehair-zigzag 
interfaees (see below), where the eleetron transits between two extreme regimes, i.e., an ultrarel- 
ativistie (i.e., ineluding the limit of vanishing earner mass) Dirae regime (armehair segment) and 
a nonrelativistie Sehrodinger regime (zigzag segment). We have thus been led to adopt the same 
Dirae-type transfer-matrix approaeh as with the armehair GNRs, but with nonvanishing potentials 
V = This amounts to shifts (in opposite senses) of the energy seales for partiele and hole 

exeitations, respeetively, and it yields the desired vanishing value for the valenee-to-eonduetion 
gap of zigzag GNRs. 

The ealeulated DFP eonduetanee that reproduees well the TB-NEGF result for the ZZZ june¬ 
tion [Fig. I^a)] is displayed in Fig. [^e); the parameters used in the DFP ealeulation are given in 
the eaption of Fig. We note that the earner mass (mi) in the eentral zigzag segment exhibits 
an energy dependenee. This is similar to a well known effeet (due to nonparabolieity in the E — k 
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dispersion) in the transport theory of usual semieonduetors.l^ We further note that the average 
mass assoeiated with a zigzag segment is an order of magnitude larger than that found for semi- 
eondueting armehair segments of similar width (see eaption in Fig.[^, and this yields energy levels 
~ {u/LY elose to the nonrelativistie limit [see Eq. (??)]. We note that the FP pattern of the ZZZ 
junetion belongs to the eategory FP-C. 

Mixed armchair-zigzag-armchair segmented GNRs. 

Fig. 1^ (right eolumn) presents an example of a mixed armehair-zigzag-armehair (AZA) junetion, 
where the eentral segment has again zigzag edge terminations [see lattiee sehematie in Fig.|^II)]. 
The eorresponding TB-NEGF eonduetanee is displayed in Fig.|^b). In spite of the different mor¬ 
phology of the edges between the leads (armehair) and the eentral segment (zigzag), the eondue¬ 
tanee profile of the AZA junetion [Fig.j^b)] is very similar to that of the ZZZ junetion [Fig.|^a)]. 
This means that the eharaeteristies of the transport are determined mainly by the eentral segment, 
with the left and right leads, whether zigzag or armehair, aeting as reservoirs supplying the im¬ 
pinging eleetrons. 

The DFP result reprodueing the TB-NEGF eonduetanee is displayed in Fig. [^d), and the pa¬ 
rameters used are given in the eaption. We stress that the mixed AZA junetion represents a rather 
unusual physieal regime, where an ultrarelativistie Dirae-Weyl massless eharge earrier (due to the 
metallie armehair GNRs in the leads) transits to a nonrelativistie massive Sehrodinger eleetron in 
the eentral segment. We note that the FP pattern of the AZA junetion belongs to the FP-C earegory. 

Aharonov-Bohm spectra of rhombic graphene rings. 

The energy of a partiele (with onedimensional momentum px) is given by the Einstein rela- 
tivistie relation E = \J[pxVpY + where ^ is the rest mass. As aforementioned, in 

armehair graphene ribbons, the mass parameter is related to the partiele-hole energy gap. A, as 
^ = A/(2v^). In relativistie quantum field theory, the mass of elementary partieles is imparted 
through interaetion with a sealar field known as the Higgs field. Aeeordingly, the mass ^ is re- 
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Figure 4: Aharonov-Bohm spectra for rhombic armchair graphene rings; the two rings that 
we consider (I and II on the left) show different atomic arrangement at the top and bottom 
corners, (a) Tight-binding spectrum for a nanoring with type-I corners and width — 12. 
(b) TB spectrum for a nanoring with type-II comers and the same width — 12. These arm¬ 
chair graphene rings are semiconducting (type-I) and metallic (type-II). The three (four) lowest- 
in-energy two-membered bands are shown. The hole states (with e < 0, not shown) are symmetric 
to the particle states (with e > 0). (c,d) DKP spectra reproducing the TB ones in (a) and (b), 
respectively. Insets in (c) and (d): schematics of the Higgs fields (position-dependent mass) ^{x) 
employed in the DKP modeling, (j) (x) is approximated by steplike functions ; i counts the three 

regions of each half of the rhombus (Lj”^ = = a and = b), and n {n= 1,2). The non-zero 

(constant) variable-mass values of ^ (x) are indicated by yellow (red) color when positive (nega- 
five). The parameters used in the DKP modeling are: (c) a = \.3ao, b = 66ao, m\' = = 0, 

= niQ = Q.l5t/vp [see corresponding schematic inset in (c)] and (d) a = 6.5ao, b = 55ao, 

= 0, — (—l)"mo with mo = 0.15t/v^ [see schematic inset in (d)]. Note the two- 

membered braided bands and the “forbidden” band [within the gap, in (b) and (d)]. The twofold 
forbidden band with e ~ 0 appears as a straight line due to the very small amplitude oscillations 
of its two members, ao = 0.246 nm is the graphene lattice constant and t = 2.7 eV is the hopping 
parameter. The edge terminations of both the inside and outside sides of the ring are armchair. 
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placed by a position-dependent Higgs field ^{x) = m{x), to whieh the relativistie fermionie field 
'P(x) eouples through the Yukawa Lagrangian^^^^^^l^^y = —(/3 being a Pauli matrix). In 
the elementary-partieles Standard Model,sueh eoupling is responsible for the masses of quarks 
and leptons. For ^{x) = (eonstant) = ^o, and the massive fermion Dirae theory is reeov- 
ered. 

We exploit the generalized Dirae physies governed by a total Lagrangian density -|- 

, where the fermionie part is given by 

+ ( 20 ) 

at ox 

and the sealar-field part has the form 

= ( 21 ) 

with the potential V (^) (seeond term) assumed to have a double-well form; ^ and 0 o are eon- 

stants. 

Heneeforth, the Dirae equation (see Methods) is generalized as 

E^' + ihvFa^-p^(x)^' = 0. ( 22 ) 

ox 

In one dimension, the fermion field is a two-eomponent spinor 'P = ( 1 / 4 , u and I stand, 
respeetively, for the upper and lower eomponent and a and /3 ean be any two of the three Pauli 
matriees. 

A graphene polygonal ring ean be viewed as made of eonneeted graphene-nanoribbon frag¬ 
ments (here we eonsider aGNRs). The exeitations of an infinite aGNR are deseribed by the ID 
massive Dirae equation, see Eq. (??) with a = O 2 , = Ci, and ^{x) = ^0 = = \h — hi- The 

two (in general) unequal hopping parameters 0 and t 2 are assoeiated with an effeetive ID tight- 
binding problem) and are given^^by 0 = — 2 teos[p;r/-I- 1 )], p = 1 , 2 ,..., and h — —t; 
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Figure 5: Wave functions for an excitation belonging to the “forbidden” solitonic band, (a) A- 

sublattice (red) and 5-sublattice (blue) components of the TB state with energy e = 0.280 x 
at = 0, belonging to the forbidden solitonic band of the type-II nanoring with = 12 [see Fig. 
I^b)] . (b) Upper (red) and lower (blue) spinor components for the corresponding state (forbidden 
band) according to the DKP spectrum [see Fig. j^d)], reproducing the TB behavior of the type- 
II nanoring with = 12 (mo = 0.15t/v^). The TB and DKP wave functions for all states 
of the solitonic band are similar to those displayed here. The wave functions here represent a 
pair of solitons. For contrast, see Fig. 10 in Ref. |2T] which portrays schematically the spinor 
for a single fermionic soliton attached to a Higgs field with a smooth kink-soliton analytic shape 
0^(x) = 0otanh (^^y^/2(j)oxJ. ^k{^) is a solution^i^I^ of the Lagrangian in Eq. (??). DKP densities 
in units of 10~^/ao. 


is the number of carbon atoms specifying the width of the nanoribbon and t = 2.7 eV is the 
hopping parameter for 2D graphene. The effective!^ TB Hamiltonian of an aGNR has a form simi¬ 
lar to that used in tranj'-polyacetylene (a single chain of carbon atoms). In tran^-polyacetylene, the 
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inequality of ti and t 2 (referred to as dimerization) is a eonsequenee of the aforementioned Peierls 
distortion indueed by the eleetron-phonon eoupling. For an armehair graphene ring, this inequality 
is a topologieal effeet assoeiated with the geometry of the edge and the width of the ribbon. We 
reeall that as a funetion of their width, , the armehair graphene nanoribbons fall into three 
elasses: (I) — 31 (semieondueting, A > 0), (II) = 3/ + 1 (semieondueting, A > 0), and 

(III) = 31 + 2 (metallie A = 0), Z = 1,2,3,.. .. 

We adapt the “erystal” approaeh^^to the Aharonov-Bohm (AB) effeet, and introduee a virtual 
Dirae-Kronig-Penneyl^ (DKP) relativistie superlattiee (see Methods). Charged fermions in a per- 
pendieular magnetie field eireulating around the ring behave like eleetrons in a spatially periodie 
strueture (period Q)) with the magnetie flux 4>/d>o (d>o = hcle) playing the role of the Bloeh wave 
veetor k, i.e., 2;rd>/4>o = [see the eosine term in Eq. (??)]. 

Naturally, nanorings with arms made of nanoribbon segments belonging to the semieondueting 
elasses may be expeeted to exhibit a partiele-hole gap (partiele-antipartiele gap in RQF theory). 
Indeed this is found for a rhombie armehair graphene ring (AGR) [see gap A in Fig. j^a)] with a 
width of = 12 earbon atoms having type-I eorners. Suprisingly, a rhombie armehair graphene 
nanoring of the same width = 12, but having eorners of type-II, demonstrates a different 
behavior, showing a “forbidden” band (with e ~ 0) in the middle of the gap region [see Fig.Qb)]. 

This behavior of rhombie armehair graphene rings with type-II eorners ean be explained through 
analogies with RQF theoretieal models, deseribing single zero-energy fermionie solitons with frae- 
tional eharge^^^l^ or their modifieations when forming soliton/anti-soliton systems (A solu¬ 
tion of the equation of motion eorresponding to Eq. (??), is a Z 2 kink soliton, The solution 

of Eq. (??) with (j) = (j)/c(x) is the fermionie soliton.) We model the rhombie ring with the use of 
a eontinuous ID Kronig-Penney^ model (see Methods) based on the generalized Dirae equation 
(??), allowing variation of the sealar field (j)(x) along the ring’s arms. We find that the DKP model 
reproduees [see Fig. |^d)] the speetrum of the type-II rhombie ring (ineluding the forbidden band) 
when eonsidering alternating masses ±mo assoeiated with eaeh half of the ring [see inset in Fig. 

Bd)]- 
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In analogy with the physics of tranj'-polyacetylene (see remarks in the introductory section), 
the positive and negative masses correspond to two degenerate domains associated with the two 
possible dimerization pattems ^^^^ ... —1\ —12 — ti —12 — ... and ... —12 — ti —12 —1\ — ..which 
are possible in a single-atom chain. The transition zones between the two domains (here two of 
the four corners of the rhombic ring) are referred to as the domain walls. 

For a single soliton, a (precise) zero-energy fermionic excitation emerges, localized at the 
domain wall. In the case of soliton-antisoliton pairs, paired energy levels with small positive and 
negative values appear within the gap. The TB spectrum in Fig.^b) exhibits a forbidden band of 
two paired -f /— levels, a property fully reproduced by the DKP model that employs two alternating 
mass domains [Fig.|^d). The twofold forbidden band with e ~ 0 appears as a straight line due to 
the very small amplitude oscillations of its two members, uq = 0.246 nm is the graphene lattice 
constant and t = 2.7 eV is the hopping parameter. 

The strong localization of a fraction of a fermion at the domain walls (two of the rhombus’ 
comers), characteristic of fermionic solitons^^and of soliton/anti-soliton pairs,l^is clearly seen in 
the TB density distributions (modulus of single-particle wave functions) displayed in Fig.j^a). The 
TB A (B) sublattice component of the tight-binding wave functions localizes at the odd numbered 
comers. These alternating localization patterns are faithfully reproduced [see Fig. |^b)] by the 
upper, t/4, and lower, t///, spinor components of the continuum DKP model. The soliton-antisoliton 
pair in Fig. [^b) generates an e/2 charge fractionization at each of the odd-numbered corners, 
which is similar to the e/2 fractionization familiar from polyacetylene. 

The absence of a forbidden band (i.e., solitonic excitations within the gap) in the spectrum of 
the type-I rhombic nanorings [see Fig. |^a)] indicates that the comers in this case do not induce 
an alternation between the two equivalent dimerized domains (represented by ±mo in the DKP 
model). Here the comers do not act as topological domain walls. Nevertheless, direct correspon¬ 
dence between the TB and DKP spectra is achieved here too by using a variable Higgs field defined 
as ^(x) = m|"^(x) with = 0 and = mg = 0.15t/v^ [see the schematic inset in Fig. 

Qc) and the DKP spectrum plotted in the same figure]. 
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Discussion 


In this paper we focused on manifestations of relativistic and/or nonrelativistic quantum behavior 
explored through theoretical considerations of transport in graphene nanostructures and spectral 
and topological effects in graphene nanorings in the presence of magnetic fields. In particular, we 
investigated the emergence of new behavior of electrons in atomically precise segmented graphene 
nanoribbons (GNRs) with different edge terminations (armchair, zigzag and mixed ones), and in 
graphene rings. To these aims we have employed tight-binding calculations of electronic states 
with and without applied magnetic fields, the non-equilibrium Green’s function transport theory, 
and a newly developed Dirac continuum model that absorbs the valence-to-conductance energy 
gaps as position-dependent masses, including topological-in-origin mass-barriers at the contacts 
between segments. 

The electronic conductance has been found to exhibit Fabry-Perot oscillations, or resonant tun¬ 
neling, associated with partial confinement and formation of a quantum box (resonant cavity) in 
the junction. Along with the familiar optical FP oscillations, exhibiting equal spacing between 
neighboring peaks, that we find for massless electrons in GNRs with metallic armchair central seg¬ 
ments, we find other FP categories that differ from the optical one. In particular, our calculations 
reveal: (a) A massive relativistic FP pattern exhibiting a valence-to-conduction gap and unequal 
peak spacings. This pattern is associated with semiconducting armchair nanoribbon central seg¬ 
ments, irrespective of whether the armchair leads are matallic or semiconducting, (b) A massive 
non-relativistic FP pattern with 1/L^ peak spacings, but with a vanishing valence-to-conduction 
gap. This pattern is the one expected for the carriers in usual semiconductors described by the 
(nonrelativistic) Schrodinger equation, and it is associated with zigzag nanoribbon central seg¬ 
ments, regardless of whether zigzag or metallic armchair leads are used. 

Perfect quantized-conductance flat steps were found only for uniform GNRs. In the absence 
of extraneous factors, like disorder, in our theoretical model, the deviations from the perfect 
quantized-conductance steps were unexpected. However, this aforementioned behavior obtained 
through TB-NEGF calculations is well accounted for by a ID contimuum fermionic Dirac-Fabry- 


29 



Perot interference theory (see Methods). This approach employs an effective position-dependent 
mass term in the Dirac Hamiltonian to absorb the finite-width (valence-to-conduction) gap in arm¬ 
chair nanoribbon segments, as well as the barriers at the interfaces between nanoribbon segments 
forming a junction. For zigzag nanoribbon segments the mass term in the Dirac equation re¬ 
flects the nonrelativistic Schrodinger-type behavior of the excitations. The carrier mass in zigzag- 
terminated GNR segments is much larger than the particle mass in semiconducting armchair- 
terminated GNR segments. Furthermore in the zigzag GNR segments (which are always char¬ 
acterized by a vanishing valence-to-conduction energy gap), the mass corresponds simply to the 
carrier mass. In the armchair GNR segments, the carrier mass endows (in addition) the segment 
with a valence-to-conduction energy gap, according to Einstein’s relativistic energy relation [see 
Eq. (1)]. 

We concluded with a brief discussion of the physics of electrons in segmented polygonal rings, 
which may be regarded as constructed by connecting GNR segments. Evaluation of the elec¬ 
tronic states in a rhombic graphene nanoring under the influence of an applied magnetic field in 
the Aharonov-Bohm regime, and their analysis with the use of a relativistic quantum-field theo¬ 
retical model, unveils development of a topological-in-origin zero-energy soliton state and charge 
fractionization. 

The above findings point to a most fundamental underlying physics, namely that the topology 
of disruptions of the regular honeycomb lattice (e.g., variable width segments, comers, edges) gen¬ 
erates a scalar-potential field (position-dependent mass, identifiedl^^^ also as a Higgs-type field), 
which when integrated into a generalized Dirac equation for the electrons provides a unifying 
framework for the analysis of transport processes through graphene segmented junctions and the 
nature of electronic states in graphene nanorings. 

With growing activities and further improvements in the areas of bottom-up fabrication and 
manipulation of atomically precisel^^^ graphene nanostmctures and the anticipated measurement 
of conductance through them, the above findings could serve as impetus and implements aiding 
the design and interpretation of future experiments. 
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